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.^ i Recently, some papers in the literature have shown that, from a bimetric 

theory of gravity, it is possible to produce massive gravitational waves which 
generate a longitudinal component in a particular polarization of the wave. After 
•y^ ' a review of previous works, in this paper the longitudinal response function 

r~^ . of interferometers for this particular polarization of the wave is computed in 

f^ ' two different gauges, showing the gauge invariance, and in its full frequency 

OO . dependence, with specific application to the Virgo and LIGO interferometers. 



1 Introduction 

The data analysis of interferometric gravitational waves (GWs) detectors has 
recently started (for the current status of GWs interferometers see [H El El lU [3 
[6l [II [8]) and the scientific community hopes in a first direct detection of GWs 
in next years. 

Detectors for GWs will be important for a better knowledge of the Universe 
and also to confirm or ruling out the physical consistency of General Relativity 
or of any other theory of gravitation [9l [TOl HH [El UHl El [15] . This is because, in 
the context of Extended Theories of Gravity, some differences between General 
Relativity and the others theories can be pointed out starting by the linearized 
theory of gravity [21 [IHl HH [HI [E] ■ In this picture, recently, some papers in the 
literature have shown that, from a bimetric theory of gravity, it is possible to 
produce massive gravitational waves which generate a longitudinal component 
in a particular polarization of the wave [14|,ll5j. After a review of previous works 
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(i.e. the work of de Paula, Miranda and Marinho [TS] and my previous research 
[M]) on this topic, which is due to provide a context to bring out the relevance of 
the results, in this paper the longitudinal response function of interferometers 
for this particular polarization of the wave is computed in its full frequency 
dependence and in two different gauges, showing the gauge invariance, with 
specific application to the Virgo and LIGO interferometers. 

2 A review of previous results on massive grav- 
itational waves from the bimetric theory of 
gravity 

An extension of linearized general relativity which takes into account massive 
gravitons gives a weak-field stress-energy tensor [T31 [H] 

T^,T^ = "if {V ~ Iii9o'r'h^0]i9o),.}, (1) 

where rUg is the mass of the graviton, and {go)p,u the non-dynamical back- 
ground metric (note: differently from [15| in this paper we work with G = 1, 
c = 1 and fi, := 1, exactly like in [H]). In this way the field equations can be 
obtained in an cinstenian form like 

G^. = -8^(T^,+Tjir^), (2) 

where T^i/ is the ordinary stress-energy tensor of the matter. General rela- 
tivity is recovered in the limit rUg -^ 0. 

Calling Qp^u the dynamic metric and putting 

gtiv = Vt^y + h^u (3) 

with |ft.^i/| <C 1 equation ([2]) can be linearized in vacuum (i.e. T^^ = 0) 
obtaining 

Uh^y = m^hf^y, (4) 

where D is the d'Alembertian operator and h^u = h^i, — -2'ntJ.v- 
The general solution of this equation is [131 [H] 

h^j.u = efiv e:>q){ik°'xa), (5) 

where e^i, is the polarization tensor. 



The condition of normalization k'^ka = rni gives k = \/^'^ ^ "i^ and a speed 
of propagation 



m^ 



vico) = ^^— , (6) 



which is exactly the velocity of a massive particle with mass mg (it is also 
the group- velocity of a wave-packet p^ [Ml [TBI [TTl fT8] ) . 

Thus, assuming that the wave is propagating in the z direction, the metric 
perturbation ([5]) can be rewritten like 

hf_iu = Gfiu cxp{ikz - iuj). (7) 

Using a tctradc formalism, the authors of [15j found six independent polar- 
izations states (see equations 28-33 of [H]), while in [T3] it has been shown that, 
from the polarization labelled with $22 in [IS] (equations 32 and 38 of [15]), a 
longitudinal force is present. 

In fact, let us consider equation 38 of [TS]. Putting hg = /loo + hss, this 
equation can be rewritten as |14| 



<^22 = ^hg{t-vz). (8) 
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Taken in to account only the <I>22 polarization in equation ([5]) one gets [14| 



hf_„^{t, z) = -hg{t - vz)r]f_,^ (9) 

and the corrispondcnt line clement is the conformally flat one |14] 

ds^ = [1 + -hg{t - vz)]{~dt'^ + dz^ + dx^ + dy^). (10) 

8 

Because the analysis on the motion of test masses is performed in a laboratory 

environment on Earth, the coordinate system in which the space-time is locally 

flat is typically used and the distance between any two points is given simply by 

the difference in their coordinates in the sense of Newtonian physics [T^ [HI [TBI 

[TTldn]- This frame is the proper reference frame of a local observer, located for 

example in the position of the beam splitter of an interferometer. In this frame 

gravitational waves manifest themself by exerting tidal forces on the masses 

(the mirror and the beam-splitter in the case of an interferometer). A detailed 

analysis of the frame of the local observer is given in ref. [19], sect. 13.6. Here 

only the more important features of this coordinate system are recalled: 

the time coordinate xq is the proper time of the observer O; 

spatial axes are centered in O; 

in the special case of zero acceleration and zero rotation the spatial coor- 
dinates Xj are the proper distances along the axes and the frame of the local 
observer reduces to a local Lorcntz frame: in this case the line element reads 

m 

ds'^ = -{dx^f + S.jdx'dx^ + 0{\x^^)dx°'dxl^. (11) 

The effect of the gravitational wave on test masses is described by the equa- 
tion 

x' = —Rq^.qX , (12) 



which is the equation for geodesic deviation in this frame. 
_ Thus, to study the effect of the massive gravitational wave on test masses, 
Rqi^q has to be computed in the proper reference frame of the local observer. 
But, because the linearized Riemann tensor R^^pa is invariant under gauge 
transformations [12l [141 [TBI [19] , it can be directly computed from eq. ([9]). 
From [12] it is: 

^pi/pCT = ■^{dp.dfjhaij + d^dah^p — dadph^^ - dp_di,hap}, (13) 

that, in the case eq. ([9|), begins 

-R070 == T77{c'"ao/ig77o7 + dodjhgSo - d"dyhgT]oo - dodohgS^}; (14) 

- dodohgS^^ = -c'j/ig for = 7. (15) 

The computation has been performed in T5]. obtaining 

pi — _ 1 )i 

-^010 — 16 "'9 

-^020 = -JE^a (16) 

The third of cqs. P^ shows that the field is not transversal. 
Infact, using eq. ([T^ it results 

y = ^hy (18) 

and 

z^-—mlhg{t-vz)z. (19) 

Then the effect of the mass is the generation of a longitudinal force (in 
addition to the transverse one). 

For a better understanding of this longitudinal force, in [14] the effect on 
test masses in the context of the geodesic deviation has been analysed. 

Following [T3j one puts 
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^0,0 = Ye ^ "^* |;i,(i-z;z) = --r,,a,'^/i, + -L,,m^/.g. (20) 

Here the transverse projector with respect to the direction of propagation of 
the GW n, defined by 
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Figure 1: a GW propagating from an arbitrary direction 



and the longitudinal projector defined by 



(21) 



L.y = niUj (22) 

have been used [Mj . In this way the geodesic deviation equation (fT2|) can be 
rewritten like 



d^ 1 



frh T- T- 



1 
TTl„ 
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^h„L, 



di2-'- i6"'-^-'^-^ lQ-9-y-^^-r (23) 

Thus it appears clear what was claimed in previous discussion: the effect 
of the mass present in the GW generates a longitudinal force proportional to 
m? which is in addition to the transverse one. But if v{uj) -^■ 1 in cq. ^ 



we get nig 



0. and the longitudinal force vanishes. Thus it is clear that the 



longitudinal mode arises from the fact that the GW does no propagate at the 
speed of light. 

In [14] it has also been analizcd the dectability of the polarization ^ com- 
puting the pattern function of a detector to this massive component. One has to 
recall that it is possible to associate to a detector a detector tensor that, for an 
interferometer with arms along the u eij directions with respect the propagating 
gravitational wave (see figure 1), is defined by [21 [Ml [17] 



2^ ^ 



(24) 



If the detector is an interferometer [Tl[2l[3l[4l[5l[6l[7l[8], the signal induced 
by a gravitational wave of a generic polarization, here labelled with s(i), is the 
phase shift, which is proportional to [21 [HI [17] 



s{t) ~ D''R,ojo (25) 

and, using equations ([20]), one gets 



s{t) sin^ 6* cos 20. (26) 

The angular dependence ()26p , is different from the two well known standard 
ones arising from general relativity which are, respectively 

(1 + cos^ 6*) cos 2(/) 
for the + polarization and 

— cos 6 sin 2-d 

for the X polarization. 

Thus, in principle, the angular dependence (|26p could be used to discriminate 
among the bimetric theory and general relativity, if present or future detectors 
will achieve a high sensitivity. 

3 The longitudinal response function 

But there is a problem. The function (j26|) is not the general form of the re- 
sponse function, but it is only a good approximation for long wavelengths (i.e. 
the wavelength of the wave is much larger than the linear dimension of the 
interferometer) [SJ [31 [TH [T7] . Now the full frequency dependent response func- 
tion will be computed. For a sake of semplicity we will consider the case of a 
massive gravitational way propagating in a direction parallel to one arm of the 
interferometer. 

We first performe the computation in the gauge (jlOp , which is not the gauge 
of the local observer in which previous computations have been performed. Eq. 
PI7)) can be rewritten as [H] 

where r is the proper time of the test masses. 

From eqs. ([T0|) and ([27|l the geodesic equations of motion for test masses 
(i.e. the beam-splitter and the mirrors of the interferometer), can be obtained 

d X r\ 



(28) 



d\ 





dh 
dT^ 


2 {l+^h,Y 


d^Z 

dr^ 


l9.(l+|'^s) 
2 (l + i/,g)2 



Note: equations ([^5]) are different from equations pT|) , ^TE\i and ([TO]) because 
the gauge pO|) is not the gauge of the local observer where equations ([17]), p8|) 
and (|19p have been performed. The gauge-invariance between the two gauges 
will be shown in Section 4. The first and the second of cqs. (|^5|) can be 
immediately integrated obtaining 

= Ci^ const. (29) 

(IT 

^^0-2= const. (30) 

UT 

In this way eq. ^7} becomes 

If we assume that test masses are at rest initially we get Ci = C2 = 0. Thus 
we see that, even if the GW arrives at test masses, we do not have motion of 
test masses within the x — y plane in this gauge. We could understand this 
directly from eq. ([TU]) because the absence of the x and of the y dependences in 
the metric implies that test masses momentum in these directions (i.e. Ci and 
C2 respectively) is conserved. This results, for example, from the fact that in 
this case the x and y coordinates do not esplicitly enter in the Hamilton-Jacobi 
equation for a test mass in a gravitational field [151 HI] • 

Now we will sec that, in presence of the GW, we have motion of test masses 
in the z direction which is the direction of the propagating wave. An analysis 
of cqs. (|28p shows that, to simplify equations, we can introduce the retarded 
and advanced time coordinates (a, 6): 

a = t — vz 

(32) 

b = t -\- vz. 



From the third and the fourth of eqs. (|^5]) we have 
d da _ db[l + \hg{a)\ 



drdT [i + lhg{a)Y 
This equation can be integrated obtaining 

da 



0. (33) 



d7="' 



(34) 

where a is an integration constant. From eqs. PT]) and (|M]) . we also get 

dh f3 



dT 1 + ^hg 



(35) 



where 3 = —, and 



T = /3a + 7, (36) 

where the integration constant 7 correspondes simply to the retarded time 
coordinate translation a. Thus, without loss of generality, we can put it equal 
to zero. Now let us see what is the meaning of the other integration constant 
f3. We can write the equation for z from eqs. ([M]) and ([55]) : 



When it is hg = (i.e. before the GW arrives at the test masses) eq. (|37p 
becomes 

But this is exactly the initial velocity of the test mass, then we have to 
choose /3 = 1 because we suppose that test masses are at rest initially. This also 
imply a = 1. 

To find the motion of a test mass in the z direction we see that from eq. 
(|36|) we have dr = da, while from eq. ((35l) we have db = _f'[, . Because it is 

vz = ^^ we obtain 



2 



dz = l(^-.a), (39) 



which can be integrated as 

_ 1 ft-vz ^hg(a) , 



(40) 



where zq is the initial position of the test mass. Now the displacement of 
the test mass in the z direction can be written as 



2v J~oo i+i/ig(a) 

- 2v J-oc l+lhg(a)""- 

We can also rewrite our results in function of the time coordinate t: 

x{t) = Xo 

y{t) = 2/0 

-(0 = ^^-htrii^/i-) 

T{t) = t-Vz[t), 



(41) 



(42) 



Calling I and L + 1 the unperturbed positions of the beam-splitter and of the 
mirror and using the third of eqs. (|42p the varying position of the beam-splitter 
and of the mirror are given by 

-"■^^ ' 2v J — oo l+-ghg(a) ^ ' 

(43) 

ZM(t)=L + /-^/!-J^^^'^4|^d(a) 

But we are interested in variations in the proper distance (time) of test 
masses, thus, in correspondence of eqs. P5)) . using the fourth of eqs. (|42p we 
get 



Then the total variation of the proper time is given by 

1 r^^"^^'^ Ih.ia) 



(44) 



A T{t) = TM(i) - TBsit) =vL-- ^^'^ d{a). (45) 

In this way, recalling that in the used units the unperturbed proper distance 
(time) is T = L, the difference between the total variation of the proper time 
in presence and the total variation of the proper time in absence of the GW is 

1 /•*-(^+') ^hg{a) 



ST{t) ^ Ar(t) L = -L{v + 1) - ;- / ' ,": . d{a). (46) 

2 Jt-vi i^+sKw 

This quantity can be computed in the frequency domain, defining the Fourier 
transform of hg as 

/CXD 
dt hg{t)exp{iujt). (47) 

and using the translation and derivation Fourier theorems, obtaining 
Sriuj) = {L(l - 1-2) exp[iwL(l + v)] + ^^^iv-.^y 

[exp[2iujL]{v + l)^{-2i + LoL{v - 1) + 2Lexp[iujL{l + v)] (48) 

{6iv + 2iv^ -ujL + ujLv^) + L{v + lf{-2i + ujL{v + l))}|/ig. 
A "signal" can be also defined: 

Siu;) ^ ^ = {(1 - v') exp[zc.L(l + vg)] + .^^^y.y. 
[eyi-p[2iujL]{v + lf{-2i + ujL{v- 1) + 2cx.p[iLuL{l + v)] (49) 

{6iv + 2iv^ -LuL + ujLv^) + L{v + \f{~2i + ujL{v + l))}\hg- 
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Figure 2: the absolute value of the longitudinal response function ((50|) of the 
Virgo interferometer (L = 3Km) to a GW arising from the bimetric theory of 
gravity and propagating with a speed of 0.1c (non relativistic case). 



Then the function 

Ti{uj) = (1 - v^) exp[iujL{l + vg)] + 2^l(J^-i)^ 

[eyip[2iuL]{v + If {-2i + ujL{v - 1) + 2exp[icjL(l + w)] (50) 

{6iv + 2m)3 -ujL + uLv^) + {v + \f{-2i + ujL{v + 1))], 

is the response function of an arm of our interferometer located in the z-axis, 
due to the longitudinal component of the massive gravitational wave arising from 
the bimetric theory of gravity and propagating in the same direction of the axis. 

For ?; ^ 1 it is Ti{uj) -> 0. 

In figures 2, 3 and 4 are shown the response functions ([5(1)) for an arm of 
the Virgo interferometer (i = 3Km) for v = 0.1 (non-relativistic case), v = 0.9 
(relativistic case) and v = 0.999 (ultra-relativistic case). We see that in the 
non-relativistic case the signal is stronger as it could be expected (for rUg — > 
we expectTi(a;) — » 0). In figures 5, 6, and 7 the same response functions are 
shown for the Ligo interferometer {L ~ AKm). 
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Figure 3: the absolute value of the longitudinal response function ([5D|) of the 
Virgo interferometer {L = 3Km) to a GW arising from the bimetric theory of 
gravity and propagating with a speed of 0.9 (relativistic case). 
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Figure 4: the absolute value of the longitudinal response function ((50|) of the 
Virgo interferometer (L = 3Km) to a GW arising from the bimetric theory of 
gravity and propagating with a speed of 0.999 (ultra relativistic case). 
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Figure 5: the absolute value of the longitudinal response function ([50]) of the 
LIGO interferometer {L = AKm) to a GW arising from the bimetric theory of 
gravity and propagating with a speed of 0.1c (non relativistic case). 
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Figure 6: the absolute value of the longitudinal response function ((50|) of the 
LIGO interferometer {L = AKm) to a GW arising from the bimetric theory of 
gravity and propagating with a speed of 0.9c (relativistic case). 
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Figure 7: the absolute value of the longitudinal response function of the LIGO 
interferometer (L = AKm) to a GW arising from the bimetric theory of gravity 
and propagating with a speed of 0.999c (ultra relativistic case). 
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4 Gauge invariance of the longitudinal response 
function 

For a sake of completeness, now the gauge invariance of the longitudinal response 
function between the gauge (fTU]) and the gauge of the local obsercr will be shown. 

Equations (fT7|) . p^ and P^ give the tidal acceleration of the test mass 
caused by the gravitational wave respectly in the x direction, in the y direction 
and in the z direction [TBI HQ] . 

Equivalently we can say that there is a gravitational potential [TH [THl [201 ■ 



ViT',t) = - — hgit^-)[x^ + y^] + —mpJ hg{t~vz)wdw, (51) 



which generates the tidal forces, and that the motion of the test mass is 
governed by the Newtonian equation 

f' = - V ^- (52) 

To obtain the longitudinal component of the gravitational wave the solution 
of eq. P^ has to be found. 

For this goal the perturbation method can be used [111 [2D] . A function of 
time for a fixed z, 'ijj{t — vz), can be defined [16j, for which it is 

Tp{t — vz) = hg{t — vz) (53) 

(note: the most general definition is Tjj{t — vz) + a{t — vz) + b, but, assuming 
only small variatons in the positions of the test masses, it results a = b = 0). 
In this way it results 

Sz{t~vz) = -—mlzoipiit-vz). (54) 

A feature of the frame of a local observer is the coordinate dependence of 
the tidal forces due by gravitational waves which can be changed with a mere 
shift of the origin of the coordinate system [16l [20] : 



X + x' , y ^ y + y' and z ^ z + z' . (55) 



The same applies to the test mass displacements, in the z direction, eq. ([54]) . 
This is an indication that the coordinates of a local observer are not simple as 
they could seem [TBI [20] ■ 

Now, let us consider the relative motion of test masses. A good way to 
analyze variations in the proper distance (time) of test masses is by means of 
"bouncing photons" (see refs. [TBI [50] and figure 8). A photon can be launched 
from the beam-splitter to be bounced back by the mirror. It will be assumed 
that both the beam-splitter and the mirror are located along the z axis of 
our coordinate system (i.e. an arm of the interferometer is in the z direction, 
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Photon 



Photon 



Figure 8: photons can be launched from the beam-spHtter to be bounced back 
by the mirror 




Figure 9: the beam sphtter and the mirror are located in the direction of the 
propagating GW 
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which is the direction of the propagating massive gravitational wave and of the 
longitudinal force, see also Figure 9). 

It will be shown that, in the frame of a local observer, two different effects 
have to be considered in the calculation of the variation of the round-trip time 
for photons, like in [1^1. Note that in [2D] the considered effects were three, but, 
if we put the beam splitter in the origin of our coordinate system, the third 
effect vanishes |16) . 

The unperturbed coordinates for the beam-splitter and the mirror are Xb = 
and Xm = L. Thus, the unperturbed propagation time between the two masses 
is 

T = L. (56) 

From cq. ([54|) it results that the displacements of the two masses under the 
influence of the gravitational wave are 

Szbit) = (57) 

and 

(58) 



(59) 



(60) 

But there is the problem that, for a large separation between the test masses 
(in the case of Virgo or LIGO the distance between the beam-splitter and the 
mirror is three or four kilometers), the definition (|59p for relative displacement 
becomes unphysical because the two test masses are taken at the same time and 
therefore cannot be in a casual connection [TBI HO] . The correct definitions for 
our bouncing photon can be written like 

SLi{t)=Sz,„it~vL)-6zbit-Ti) (61) 

and 

SL2it) = Sz„,{t -vL- T2) - 5zb{t), (62) 

where Ti and T2 are the photon propagation times for the forward and return 
trip correspondingly. According to the new definitions, the displacement of one 
test mass is compared with the displacement of the other at a later time to 
allow for finite delay from the light propagation. Note that the propagation 
times Ti and T2 in eqs. (j6ip and ([62]) can be replaced with the nominal value T 
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6z„T,{t - vL) ^ -~—m'^gLtlj{t - vL) 


In this way, the 


relative displacement, is 


5L{t) 


= 5zyn{t - vL) - 5zb{t) = -—rn^gL^I; 


Thus it results 






5L{t) 5T{t) 1 2,,, r^ 
L - T - 16™^^(* ^'^' 



because the test mass displacements are alredy first order in hg. Thus, for the 
total change in the distance between the beam splitter and the mirror in one 
round-trip of the photon, it is 



6Lr.tXt) =SLi{t-T) + SL2{t) = 26z.,n{t-vL-T)-Szb{t)-Szb{t-2T), (63) 
and in terms of ip and of the mass of the gravitational wave: 

6Lr.t.{t) = -^mlL,Pit-vL-T). (64) 

The change in distance (|64p leads to changes in the round-trip time for 
photons propagating between the beam-splitter and the mirror: 

^ = -im>(t-.L-r). (65) 

In the last calculation (variations in the photon round-trip time which come 
from the motion of the test masses inducted by the massive gravitational wave), 
it was implicitly assumed that the propagation of the photon between the beam- 
splitter and the mirror of our interferometer is uniform as if it were moving in a 
flat space-time. But the presence of the tidal forces indicates that the space-time 
is curved. As a result another effect after the previous has to be considered, 
which requires spacial separation [HI [20] . 

For this effect we consider the interval for photons propagating along the 
z-axis 

ds^ ^ goodt^ + dz"^ . (66) 

The condition for a null trajectory {ds — 0) gives the coordinate velocity of 
the photons 

«^^(g)2 = l + 2y(f,z), (67) 

which to first order in hg is approximated by 

Vf^±[\ + V{t,z)], (68) 

with + and — for the forward and return trip respectively. Knowing the 
coordinate velocity of the photon, the propagation time for its travelling between 
the beam-splitter and the mirror can be defined: 

Ti{t) = / - (69) 



and 



ut) = r^'^ ^^. (70) 
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The calculations of these integrals would be complicated because the bound- 
ary z„i{t) is changing with time. In fact it is 

zb(i) = 6zb{t) = (71) 

but 

Z^{t)=L + 5Zrn{t). (72) 

But, to first order in /ig, this contribution can be approximated by SL2{t) 
(see eq. (|5^ ). Thus, the combined effect of the varying boundary is given 
by SiT(t) in eq. (j65p . Then only the times for photon propagation between 
the fixed boundaries and L have to be calculated. Such propagation times 
will be denoted with ATi,2 to distinguish from Ti.2. In the forward trip, the 
propagation time between the fixed limits is 

dz 



where t' is the retardation time which corresponds to the unperturbed pho- 
ton trajectory: 

t' = t-{L-z) 

(i.e. t is the time at which the photon arrives in the position i, so L — z = 
t-t'). 

Similiary, the propagation time in the return trip is 

AT2(i)=T- / V{t\z)dz, (74) 

where now the retardation time is given by 

t' ^t- z. 

The sum of ATi{t — T) and /S.T2{t) gives the round-trip time for photons 
traveling between the fixed boundaries. Then the deviation of this round-trip 
time (distance) from its unperturbed value 2T is 

62T{t)= f [V{t^2T + z,z) + V{t-z,z)]dz. (75) 

Jo 



From eqs. (|5T|) and ([75]) it results: 

52T{t) ~ i?77g Jq [/q \hg{t — 2T + w — vw)wdw + J^ \hg{t — w — vw)wdw]dz - 

= \ral \^\\hg(t ^vz-2T + z) + \hg{t - vz - z)]z^dz+ 

— jtri^g Jp [Jg \h'g{t — 2T + w ~ vw)z'^dw + Jg \h'g{t — w — v'w)z'^d'w\dz, 

(76) 
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Thus the total round-trip proper distance in presence of the massive gravi- 
tational wave is: 

T = 2T + 5iT + S2T. (77) 

Now, to obtain the interferometer response function of the massive gravita- 
tional wave, the analysis will be transled in the frequency domine. 
Using the Fourier transform of ip defined from 



ip{uj)^ dt'ip{t)exp{iujt), (78) 

J —00 

eq. (j65p can be rewritten like: 

'41^ = ^ImlTU.m.) (79) 

with 

TI{lu) = exp[iw(l + vg)L]. (80) 

But. from a theorem about Fourier transforms, it is simple to obtain: 

^M = -H^, (81) 

where the Fourier transform of hg is given by equation (j47p . 
Then it results: 

Sifjuj) ml ~ 

J, == ■^^iHhgK^)^ (82) 

and, defining: 



we obtain: 



Ti^^T*(c.)-(l-«^)TJ(c.), (83) 



'-^ = \T,{.)hg[.). (84) 

On the other hand eq. (|76p can be rewritten in the frequency space like: 

Jaf (cj) = ,^(,.^1). [eM'^iuL\{v + lf{-2i + ujL{v - 1) + 

+ 2eyij>[iuuL{l + v)\{Qiv + 2iv'^ -ujL + u}Lv^)+ (85) 

+ {v + \f{-2i + ujL{v + l))]^. 
Now 
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52f{uj) _ ^ ,,,^hg{uj) 



T 
can be put, with 



T2(c.)^^, (86) 



T^2(^) ^ 2ujL(v'^-ip [exp[2zcj£](t. + l)3(-2z + ujL{v - 1) + 

2exp[iLjL(l + z;)](6w + 2iv^ - uL + ujLv'^)+ (87) 

+ (^. + l)3(-2^ + wi(^. + l))]. 
Because it is 

TKc^) = Ti(cc>) + T2(c^), (88) 

from eqs. ([80]) . (|83|) and (|87p it results that the function 

T;(w) = (1 - u^) exp[iwL(l + vg)] + 2uL(l-^-iy^ 
[exp[2iwL](i; + lY{-2i + ujL{v - 1) + 2exp[icjL(l + w)] 



(6w + 2m)3 - wL + tjLv-*) + (u + lf{-2i + ujL{v + 1))], 

is the longitudinal response function of an arm of the interferometer located 
in the z-axis, due to the longitudinal component of the massive gravitational 
wave propagating in the same direction of the axis, and one can see that equation 
([M|) is equal to equation ((501) . 

Thus, we have shown that the longitudinal response function of an arm of 
an interferometer located in the z-axis is the same in both the local Lorentz 
gauge and in the gauge PU)) . 

5 Conclusions 

This paper is an integration of previous research on massive gravitational waves 
from a bimetric theory of gravity. In the literature about this issue, it has been 
shown that massive gravitational waves arising from such a bimetric theory 
can generate a longitudinal component in a particular polarization of the wave 
[141 I15| . After a review of previous works, which was due for completeness 
and for a better understanding of the analysis, in this paper the longitudinal 
response function of interferometers for this particular polarization of the wave 
has been computed in two different gauges, showing the gauge invariance, and in 
its full frequency dependence, with specific application to the Virgo and LIGO 
interferometers . 
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